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Abstract 

The combination of large per-photon optical force and small motional mass attainable in nanocavity op- 
tomechanical systems results in strong dynamical back-action between mechanical motion and the cavity 
light field. In this work we study the optical control of mechanical motion within two different nanocavity 
structures, a zipper nanobeam photonic crystal cavity and a double-microdisk whispering-gallery resonator. 
The strong optical gradient force within these cavities is shown to introduce signifcant optical rigidity 
into the structure, with the dressed mechanical states renormalized into optically-bright and optically-dark 
modes of motion. With the addition of internal mechanical coupling between mechanical modes, a form 
of optically-controlled mechanical transparency is demonstrated in analogy to electromagnetic ally induced 
transparency of three-level atomic media. Based upon these measurements, a proposal for coherently trans- 
ferring RF/microwave signals between the optical field and a long-lived dark mechanical state is described. 

PACS numbers: 



*Electronic address: 



opainter@caltech.edu 



URL: 



http : //copilot . caltech . edu 



The coherent mixing of multiple excitation pathways provides the underlying mechanism for 
many physical phenomena. Well-known examples include the Fano resonance [1J and electro- 
magnetically induced transparency (EIT) [2], arising from the interference between excitations of 
discrete states and/or a continuum background. In the past few decades, Fano-like or EIT-like 
resonances have been discovered in a variety of physical systems, such as electron transport in 
quantum wells/dots [3, 4], phonon interactions in solids (HO, inversion-free lasers EUl, coupled 
photonic microcavities &\ QUI El [IH, and plasmonic metamaterials lfT3l . Here we report a new 
class of coherent excitation mixing which appears in the mechanical degree of freedom of nano- 
optomechanical systems (NOMS). We use two canonical systems, coupled microdisks and cou- 
pled photonic-crystal nanobeams, to show that the large optical stiffening introduced by the optical 
gradient force actuates significant coherent mixing of mechanical excitations, not only leading to 
renormalization of the mechanical modes, but also producing Fano-like and EIT-like optomechan- 
ical interference, both of which are fully tunable by optical means. The demonstrated phenomena 
introduce the possibility for classical/quantum information processing via optomechanical sys- 
tems, providing an on-chip platform for tunable optical buffering, storage, and photonic-phononic 
quantum state transfer. 

Light forces within micro-mechanical systems have attracted considerable interest of late 
due to the demonstration of all-optical amplification and self-cooling of mesoscopic mechanical 
resonators [fT4l [T5l [T6l [P71 [T8l . This technique for sensing and control of mechanical motion relies 
on the radiation pressure forces that build up in a mechanically compliant, high-Finesse optical 
cavity, resulting in strong dynamical back-action between the cavity field and mechanical mo- 
tion. More recently [|T9l l20l I2T1 l22l [23l |24j, it has been realized that guided wave nanostructures 
can also be used to generate extremely large per-photon optical forces via the gradient optical 
force [|25l . The combination of tailorable mechanical geometry, small motional mass, and large 
per-photon force in such nanostructures results in a regime of operation in which the dynamic 
response of the coupled optomechanical system can significantly differ from that of the bare me- 
chanical structure. In particular, the mechanical motion can be renormalized by the optical spring 
effect|22l |26[ |27l |28l |29l |30), creating a highly anistropic, intensity-dependent effective elastic 
modulus of the optomechanical structure. 

In this work we focus on two specific implementations of nanoscale cavity optomechanical 
systems, shown in Fig.[T] in which dynamical back-action effects are particularly strong. The first 
system consists of two patterned nanobeams in the near-field of each other, forming what has been 
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FIG. 1: a, Schematic and b, zoomed-in scanning electron microscopic (SEM) image of the double-disk 
NOMS. c, FEM-simulated electric field intensity of a transverse-electric (TE) polarized, bonded (even par- 
ity) whispering-gallery supermode between the two microdisks (shown in cross-section and for resonance 
wavelength X c 1550 nm). The double-disk bonded supermode has an optomechanical coupling coefficient 
of goivi/27r ~ 33 GHz/nm. The device studied here has a measured resonance wavelength of X c = 1538 nm 
and an intrinsic and loaded quality (Q) factor of 1.07 x 10 6 and 0.7 x 10 6 , respectively, d, Schematic, 
e, SEM image, and f, FEM-simulated bonded (even parity) optical supermode of the zipper cavity. The 
zipper cavity bonded supermode has an optomechanical coupling coefficient of goyi/ln 68 GHz/nm, a 
measured resonance wavelength of X c = 1545 nm, and an intrinsic and loaded Q-factor of 3.0 x 10 4 and 
2.8 x 10 4 , respectively. Additional details for both devices are in the Appendices and in Refs. I22ll23l . 

termed a zipper cavity Il22l [3TTl . In this cavity structure the patterning of the nanobeams localizes 
light through Bragg-scattering, resulting in a series of high Finesse (J? ^3x 10 4 ), near-infrared 
(A ~ 1550 nm) optical supermodes of the beam pair. Clamping to the substrate at either end of the 
suspended beams results in a fundamental in-plane mechanical beam resonance of frequency ~ 8 
MHz. The second cavity optomechanical system is based upon a whispering-gallery microdisk 
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optical cavity structure. By creating a pair of microdisks, one on top of the other with a nanoscale 
gap in between, strong optical gradient forces may be generated between the microdisks while 
maintaining the benefits of the low-loss, high-<2 (Q > 10 6 ) character of the whispering-gallery 
cavity. As shown schematically in Fig. [T^, the double-disk structure[23J is supported and pinned 
at its center, allowing the perimeter of the disks to vibrate in myriad of different ways. Of particular 
interest in both the zipper and double-disk systems is the differential motion of the nanobeams or 
disks, in which the modulation of the gap between the elements creates a large dispersive shift in 
the internally propagating cavity light field. It is this type of motion that is strongly coupled to the 
light field, and for which the dynamical back-action is strongest. 

We begin with an analysis of the zipper cavity, in which the strong optically-induced rigidity as- 
sociated with differential in-plane motion of the nanobeams results in a dressing of the mechanical 
motion by the light field. As described in the App.|B]and in Ref. [22j, an optical fiber nanoprobe is 
couple light into and out of the zipper cavity. Optical excitation provides both a means to transduce 
mechanical motion (which is imparted on the transmitted light field through phase and intensity 
modulation) and to apply an optical-intensity-dependent mechanical rigidity via the strong optical 
gradient force. By fitting a Lorentzian to the two lowest-order in-plane mechanical resonances in 
the radio-frequency (RF) optical transmission spectrum, we display in Fig. [2^ and b the resonance 
frequency and resonance linewidth, respectively, of the two coupled mechanical modes of the 
nanobeam pair as a function of laser-cavity detuning. At large detuning (low intra-cavity photon 
number) the nanobeams' motion is transduced without inducing significant optical rigidity, and the 
measured mechanical resonances are split by ~ 200 kHz, with similar linewidths (damping) and 
transduced amplitudes (Fig. [2jc)). As the laser is tuned into resonance from the blue-side of the 
cavity, and the intra-cavity photon number increases (to ~ 7000), the higher frequency resonance 
is seen to significantly increase in frequency while the lower frequency mode tunes to the aver- 
age of the independent beam frequencies with its transduced amplitude significantly weaker. The 
linewidth of the high frequency resonance also tends to increase, while that of the lower frequency 
mode drops. Tuning from the red-side of the cavity resonance reverses the sign of the frequency 
shifts and the roles of the high and low frequency modes. 

A qualitative understanding of the light-induced tuning and damping of the zipper cavity 
nanobeam motion emerges if one considers the effects of squeeze-film damping [|32l . Squeeze- 
film effects, a result of trapped gas in-between the beams (measurements were performed in 1 atm. 
of nitrogen), tend to strongly dampen differential motion of the beams and should be negligible 
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FIG. 2: a, Mechanical frequency and b, linewidth of the fundamental in-plane mechanical resonances of 
the zipper cavity's coupled nanobeams as a function of laser frequency detuning. The input power for these 
measurements is 127 /iW, corresponding to a maximum cavity photon number of ~ 7000 on resonance. 
The circles show the experimental data and the solid curves correspond to a fit to the data using Eq. fl5] ). 
Optically-transduced RF spectrum at a laser-cavity detuning of c, Ao/T f = 2.1 and d, Ao/T t = 0.32. The 
two nanobeams vibrate independently when the laser-cavity detuning is large, but are renormalized to the 
cooperative e, differential and f, common motions near resonance. 

for common motion of the beams. Similarly, the optical gradient force acts most strongly on the 
differential beam motion and negligibly on the common-mode motion. The sign of the resulting 
optical spring is positive for blue detuning and negative for red detuning from the cavity reso- 
nance. Putting all of this together, a consistent picture emerges from the data in Fig.[2]in which the 
nanobeams start out at large detuning moving independently with similar damping (the frequency 
splitting of ~ 200 KHz is attributable to fabrication assymetries in the beams). As the detuning is 
reduced, and approaches the cavity half-linewidth, the motion of the nanobeams is dressed by the 
internal cavity field into differential motion with a large additional optical spring constant (either 
positive or negative) and large squeeze-film damping component, and common motion with re- 
duced squeeze-film damping and minimal coupling to the light field. Due to the strong light-field 
coupling of the differential mode and the correspondingly weak coupling of the common mode, 
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we term these dressed motional states optically -bright and optically-dark, respectively. 

A quantitative model of the dressed system can be obtained by considering the following set of 
coupled equations for the mechanical motion of the nanobeams: 

d 2 xi dxi 2 F th +F +F q 

- w+ r m — + a mlXl = , (i) 

dxo ^ dxo F t h—F —F a 

^ + r m ^ + al 2 ^—-^, <2) 

Here F a is the optical gradient force, F q is the viscous force from the squeeze-film effect, and F t ^ is 
the Langevin force due to coupling of the beams to a thermal bath (further details of the properties 
of these three forces are described in the Appendices). The parameters xj and Q, m j (j = 1,2) are 
the mechanical displacement and the resonance frequency of the uncoupled, individual beams. 
We have assumed, for simplicity, that the bare mass (m), intrinsic mechanical damping (T ra ), and 
fluctuating thermal force (F t h) are the same for each of the nanobeams. The resulting spectral 
intensity of the thermally-excited optically-bright (xt, = x\— X2) differential beam motion can be 
shown to be (App.|I]), 

5^) = *^ l^)| 2 + l^)| 2 „ (3) 

m b |Li(n)Za(Q)-£ [f (Q)/m b + iT q n] [Li(fi) +L 2 (Q)) \ 

where Lj(Q.) = QjL- — Q. 2 — iT m Q. (j = 1,2), mt, = m/2 10.75 pg) is the effective mass of the 
differential mode, T is the bath temperature, and T q represents the damping rate introduced by 
the squeezed gas film in between the beams. f (Q) represents the modification to the mechanical 
susceptibility of the differential beam motion by the optical gradient force, and is approximately 
given in the bad-cavity limit (T t Q. m ) by, 

Ma) « - (tf-W 1 ) ( (Ao ^ r ^ 2)2 ) , (4) 

where goM (= dco c /dxb ~ 2^(68 GHz/nm)) is the optomechanical coupling coefficient, |ao| 2 is 
the time-averaged optical cavity energy, and T t 2^(6.9 GHz)) is the energy decay rate of the 
loaded optical cavity resonance. Lorentzian fits to the double resonances of S Xb (£l) yields the 
frequency and damping curves shown as solid lines in Fig. [2J3, displaying excellent agreement 
with the experimental measurements. 

A similar optically-induced renormalization mechanism applies to the double-disk cavity struc- 
ture shown in Fig. [T^-c. In this case, the large optical spring effect for the differential motion of 
the two microdisks excites another, more intriguing form of coherent optomechanical mixing with 



the optically dark common mode of the disks. Unlike in the zipper cavity, FEM modeling of the 
mechanics of the double-disk structure indicates a significant frequency splitting between the dif- 
ferential and common modes of motion of the double disk (shown in Fig. [3}) and c), primarily due 
to the difference in the extent of the undercut between the disk layers and the extent of the central 
pedestal which pins the two disk layers. The result is that the differential, or "flapping" motion, 
of the undercut disk region has a lower frequency of 7.95 MHz, whereas the common motion of 
the disks results in a higher frequency (14.2 MHz) "breathing" motion of the entire double-disk 
structure. 

The RF-spectrum of the transmitted optical intensity through a double-disk cavity, measured 
using the same fiber probing technique as for the zipper cavity (App. [B]), is shown in Fig. |2j: 
versus laser-cavity detuning. For the largest detuning (in which the optical spring is negligible) 
the spectrum shows a broad (2.1 MHz) resonance at 8.3 MHz and a much narrower (0.11 MHz) 
resonance at 13.6 MHz, in good corresponce with the expected frequencies of the flapping and 
breathing modes, respectively. The difference in damping between the two resonances can be 
attributed to the strong squeeze-film damping of the differential flapping motion of the disks. As 
shown in Fig. [2];, the flapping mode can be tuned in frequency via the optical spring effect from 
its bare value of 8.3 MHz all the way out to 15.7 MHz (optical input power of = 315 juW). 
In the process, the flapping mode is tuned across the breathing mode at 13.6 MHz. Although 
the optically-dark breathing mode is barely visible in the tranduced spectrum at large laser-cavity 
detunings, its spectral amplitude is considerably enhanced as the optically-bright flapping mode 
is tuned into resonance. In addition, a strong Fano-like lineshape, with ~13 dB anti-resonance, 
appears in the power spectrum near resonance of the two modes (Fig.[3]f, g, and h). 

As shown schematically in Fig. |4ja), the Fano-like interference in the optically -bright power 
spectral density can be attributed to an internal mechanical coupling between the flapping and 
breathing mechanical modes. This is quite similar to the phonon-phonon interaction during the 
structural phase transition in solids [51 [331 [341 EH [361 [37]|, m which the internal coupling between 
phonon modes produces Fano-like resonances in the Raman-scattering spectra. A Hamiltonian for 
the coupled optomechanical system can be written as, 

'^ m = in + \ mh£ ^> + it d + \ mdS **«& + KXbXd ^ (5) 

where Xj, pj, Q.j, and m 7 (j = b, d) are the mechanical displacement, kinetic momentum, intrinsic 
mechanical frequency, and effective motional mass, respectively, for the optically-bright flapping 
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FIG. 3: a,b, FEM simulated mechanical motion of the differential flapping mode (a) and the common 
breathing mode (b), with simulated frequencies of 7.95 and 14.2 MHz. The color map indicates the relative 
magnitude (exaggerated) of the mechanical displacement, c, Recorded power spectral density (PSD) of 
the cavity transmission for the double-disk, with an input power of 315 |iW. Each curve corresponds to 
a normalized laser-cavity frequency detuning, An/T, indicated in e. For display purposes, each curve is 
relatively shifted by 10 dB in the vertical axis, d, The corresponding theoretical PSD. f-h, Detailed PSD at 
three frequency detunings indicated by the arrows in e, with the experimental and theoretical spectra in blue 
and red, respectively. 



(j = b) and optically-dark breathing (j = d) mechanical modes. K represents the internal mechan- 



ical coupling between the two modes. Equation (El I together with the three external actuation 
forces (optical, squeeze-film, and thermal in nature) leads to the following equations of motion for 
the two mechanical modes: 



d 2 xt, 



dxb 



K 



,. 2 + T mb-7- + &mb x b H *d = 1 • 

dt L dt ni}, m b m b 

d 2 Xd „ d%d 2 K Fd 

dt z dt m d m d 



(6) 
(7) 



where Fj (j = b, d) represents the thermal Langevin forces and we have incorporated the squeeze- 



film damping into T m b for the optically-bright differential motion. Equations (E2) and (E3 ) lead 
to a spectral intensity for %b of (App.[E]): 

2k B Tri 4 r md + r mb \L d (Q) \ 2 
m \L b (a)L d (Q.)-rj 4 \ 



b Xb {il) - — ^ w — (8) 



where T] = K 2 /(m[,m d ) represents the mechanical coupling coefficient, Lb(Gl) = Q? — Q?, — 
iT m i ? Q.—f (Q.)/mb, and L d (Q.) =Q? — — iT md £l. As discussed in more detail in Ref. [23], the 
parameters for the double-disk are nib = 264 pg, go^/ln = 33 GHz/nm, and r f /27r = 280 MHz. 
The power spectral density of the cavity transmission, which we measure optically, is linearly 
proportional to S Xb (Q.). As shown in Fig.|2p, the theoretical model (with fitting parameter r]/2it = 
3.32 MHz) provides an excellent description of the observed phenomena. In App. |H| we present 
similar coherent mode-mixing involving two optically-dark modes and a single optically-bright 
mode in the zipper cavity. 



The mechanical response given by Eq. (E14) is directly analogous to the atomic response in 
EIT [38J. Just as in EIT, one can understand the resulting Fano lineshape in two different ways. 
The first perspective considers the interference associated with multiple excitation pathways. In 
the optomechanical system, the mechanical motion of the flapping mode is thermally excited along 
two different pathways, either directly into the broadband (lossy) flapping mode, or indirectly, 
through the flapping mode, into the long-lived breathing mode, and then back again into the flap- 
ping mode. The two excitation pathways interfere with each other, resulting in the Fano-like reso- 
nance in the spectral response of the optically bright flapping mode. An alternative, but perfectly 
equivalent view of the coupled optomechanical system considers the dressed states resulting from 
the internal mechanical coupling. In this picture the internal mechanical coupling renormalizes the 
broadband flapping mode and the narrowband breathing mode into two dressed mechanical modes, 
both broadband and optically-bright. In particular, when the flapping and breathing mechanical 
frequencies coincide, the two dressed modes are excited with equal amplitude and opposite phase 
at the center frequency between the split dressed states. Destructive interference results, suppress- 
ing excitation of the mechanical system at the line center. Consequently, the mechanical motion 
becomes purely a trapped mechanically-dark state, transparent to external excitation. As shown 
in Fig. [3^, this induced mechanical transparency is a direct analogue to EIT in atomic systems 
|2l[38l[39llini, in which the quantum interference between the transition pathways to the dressed 
states of the excited electronic state, through either |1) <-> |+) or |1) <-> |— ), leads to an induced 
spectral window of optical transparency. 
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FIG. 4: a, Schematic of an equivalent Fabry-Perot cavity system showing mechanical mode mixing. The 
mechanical motion of the cavity mirror (mi, equivalent to the optically-bright flapping mode) is primarily 
actuated by the spring k\ and the optical force. It is internally coupled to a second mass-spring system (ni2, 
equivalent to the breathing mode) actuated by the spring &2 which is decoupled from the optical wave. The 
two masses are internally coupled via spring k\2- b, A photonic analogue to the optomechanical system in- 
volving coupled resonators. Microcavity 1 is directly coupled to the external optical waveguide (equivalent 
to the optically-bright flapping mode) and also internally coupled the narrowband cavity 2 (equivalent to 
the optically-dark breathing mode), c, State diagram of an EIT-like medium. The excited state (|2)) is split 
by the optical control beam into two broadband dressed states (|+) and |— )). The dipole transition between 
ground-states 1 1) and |3) is forbidden, d, The state diagram corresponding to the optomechanical system of 
a, where |1) is the phonon vacuum state, and |2) and |3) correspond to the flapping and breathing modes, 
respectively. 

Despite the intriguing similarities between the optomechanical system studied here and EIT in 
atomic media, there are some important, subtle differences. For instance, in the optomechanical 
system, rather than the linear dipole transition of EIT, the interaction corresponds to a second-order 
transition. The dynamic backaction between the cavity field and mechanical motion creates Stokes 
and anti-Stokes optical sidebands, whose beating with the fundamental optical wave resonates 
with the mechanical motion to create/annihilate phonons (see Fig.[3]l). Functionally, this is like 
coherent Stokes and anti-Stokes Raman scattering, albeit with unbalanced scattering amplitudes 
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resulting from the coloring of the electromagnetic density of states by the optical cavity (App.|G]). 
Therefore, in analogy to EIT, it is the modulation signal carried by the incident optical wave (radio- 
frequency or microwave photons) that fundamentally probes/excites the mechanical motion and to 
which the trapped mechanically-dark state becomes transparent. Moreover, rather than tuning 
the Rabi- splitting through the intensity of a control beam resonant with the |3) <->• |2) electronic 
transition (Fig. [3]:), this optically-induced mechanical transparency is controlled via optical spring 
tuning of the resonance frequency of the optically bright flapping mechanical mode. Perhaps the 
most apt analogy to the optomechanical system can be made to the photonic resonator system 
shown in Fig. [3}i. The interference in this case is between the two optical pathways composed of 
the waveguide-coupled low-<2 optical resonator 1, and the waveguide-decoupled high-<2 resonator 
2. This interference again leads to a Fano-like resonance, or what has been termed coupled- 
resonator-induced transparency, in the optical cavity transmission [|9l [T0l[TTl[T2ll . 

Although the studies considered here involve thermal excitation of the optomechanical system, 
the same phenomena can be excited more efficiently, and with greater control, using external opti- 
cal means (App.|F]). As such, beyond the interesting physics of these devices, exciting application 
in RF/microwave photonics and quantum optomechanics exist. Similar to the information storage 
realized through EIT [|38ll4Tll42l . optical information can be stored and buffered in the dark me- 
chanical degree of freedom in the demonstrated NOMS. This can be realized through a procedure 
similar to that recently proposed for coupled optical resonators ll43l 1441 in which dynamic, adia- 
batic tuning of optical resonances are used to slow, store, and retrieve optical pulses. The corre- 
sponding optomechanical system would consist of an array of double-disk resonators, all coupled 
to a common optical bus waveguide into which an optical signal carrying RF/microwave informa- 
tion would be launched. In this scheme, a second control optical beam would adiabatically tune the 
frequency of the optically-bright flapping mode of each resonator, allowing for the RF/microwave 
signal to be coherently stored in (released from) the long-lived breathing mode through adiabatic 
compression (expansion) of the mechanical bandwidth [|43ll44ll . In comparison to the all-photonic 
system, optomechanical systems have several advantages, primarily related to the attainable life- 
time of the dark mechanical state. For example, the radial breathing mechanical mode of a similar 
whispering-gallery cavity has been shown to exhibit a lifetime of more than 2 ms ll45ll . a timescale 
more than seven orders of magnitude longer than that in demonstrated photonic-coupled-resonator 
systems [46] and comparable with EIT media ll4Tl l42l . Moreover, mechanical lifetimes of more 
than one second have recently been demonstrated using stressed silicon nitride nanobeam [47] and 
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nanomembrane [48J mechanical resonators operating in the MHz frequency regime. In the quan- 
tum realm, such a system operating in the good-cavity or sideband-resolved regime (by increasing 
either the optical Q factor [|49l or the mechanical frequency), would reduce the simultaneous cre- 
ation and annihilation of Stokes and anti-Stokes photons, enabling efficient information storage 
and retrieval at the single-quanta-level suitable for quantum state transfer. 
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APPENDIX A: DEVICE FABRICATION 

The zipper cavity is formed from a thin-film (400 nm) of tensile-stressed, stochiometric Si3N4 
deposited by low-pressure chemical vapor deposition on a silicon substrate. Electron beam- 
lithography, followed by a series of plasma and wet chemical etches, are used to form the released 
nanobeam structure. The double-disk structure is formed from a 158 nm sacrifical amorphous sil- 
icon layer sandwiched in between two 340 nm thick silica glass layers, all of which are deposited 
via plasma-enhanced chemical vapor deposition. A high temperature (1050 K) thermal anneal is 
used to improve the optical quality of the as-deposited silica layers. The microdisk pattern was 
fabricated by reactive ion etching, and the sandwiched a-Si layer was undercut by 6 /im from 
the disk edge using a sulfur hexafluoride dry release etch. This etch simultaneously undercuts the 
silicon substrate to form the underlying silicon pedestal. The final air-gap between the silica disks 
size is measured to be 138 nm due to shrinkage of the amorphous silicon layer during annealing. 

APPENDIX B: OPTICAL PROBING AND RF SPECTRUM MEASUREMENTS 

A fiber-taper optical coupling technique is used to in-couple and out-couple light from the 
zipper and double-disk cavities. The fiber taper, with extremely low-loss (88% transmission ef- 
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ficiency), is put in contact with the substrate near the cavities in order to mechanically anchor it 
during all measurements (thus avoiding power-dependent movement of the taper due to thermal 
and/or optical forces). An optical fiber polarization controller, consisting of a series of circular 
loops of fiber, is used to selectively excite the transverse-electric polarized optical modes of both 
cavities. 

RF spectra are measured by direct detection of the optical power transmitted through the cavi- 
ties using a 125 MHz bandwidth photoreceiver (noise-equivalent-power NEP= 2.5 pW/Hz 1 / 2 from 
0-10 MHz and 22.5 pW/Hz 1 / 2 from 10-200 MHz, responsivity R = 1 AAV, transimpedance gain 
G = 4 x 10 4 V/A) and a high-speed oscilloscope (2 Gs/s sampling rate and 1 GHz bandwidth). A 
pair of "dueling" calibrated optical attenuators are used before and after the cavities in order to 
vary the input power to the cavity while keeping the detected optical power level constant. The 
measured electrical noise floor is set by the circuit noise of the photodetector for the optical power 
levels considered in this work, corresponding to —125 dBm/Hz near 10 MHz. 



APPENDIX C: INTRACAVITY FIELD IN THE PRESENCE OF OPTOMECHANICAL COU- 
PLING 



In the presence of optomechanical coupling, the optical field inside the cavity satisfies the 
following equation: 

— = O'Ao - T t /2 - ig om x b )a + iy/TeAm, (CI) 

where a is the optical field of the cavity mode, normalized such that U = \a\ 2 represents the mode 
energy, and A,„ is the input optical wave, normalized such that P\ n = |A ;n | 2 represents the input 
power. r r is the photon decay rate for the loaded cavity and Y e is the photon escape rate associated 
with the external coupling. Ao = CO — (Oq is the frequency detuning from the input wave to the 
cavity resonance. g om is the optomechanical coupling coefficient associated with the optically 



bright mode, with a mechanical displacement given by x\,. In Eq. (CI ), we have neglected the 
optomechanical coupling to the optically dark mode because of its negligible magnitude. 

Well below the threshold of mechanical oscillation, the mechanical motion is generally small, 
and its impact on the intracavity optical field can be treated as a small perturbation. As a result, the 
intracavity field can be written as a(t) ps ao(t) + 8a(t), where ao is the cavity field in the absence 
of optomechanical coupling and da is the perturbation induced by the mechanical motion. They 
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satisfy the following two equations: 

dag 

dt 
dda 



dt 



(iAo - T t /2)a + i\ff e A in , 
(z'A - r,/ 2) 8a - ig om x b a . 



In the case of a continuous-wave input, Eq. (jC2J) leads to a steady state given by 

and Eq. ( |C3| ) provides a spectral response for the perturbed field amplitude of 

s~m\ ig om aox b (Q.) 
oa(Ll) - 



(C2) 
(C3) 



(C4) 



(C5) 



i(Ao + a)-^ t /2 , 

where 8a(Q) is the Fourier transform of 8a(t) defined as 8a(Q) = /+J 8a{t)e iSlt dt. Similarly, 
Xb(Q.) is the Fourier transform of x b (t). 



APPENDIX D: THE POWER SPECTRAL DENSITY OF THE CAVITY TRANSMISSION 

From the discussion in the previous section, the transmitted optical power from the cavity is 
given by 



Pt 



|A | 2 + iy/Te (AqSci -A 8a* 



(Dl) 



where Ao is the steady-state cavity transmission in the absence of optomechanical coupling. It is 
given by 



A = A t , 



(r -r e )/2-iAo 



(D2) 



r,/2-iAo ' 

where To is the photon decay rate of the intrinsic cavity. It is easy to show that the averaged cavity 
transmission is given by (Pt) = |Aq| 2 , as expected. By using Eqs. (C5), (Dl ), and (D2), we find 



the power fluctuations, 8Pr(t) = Pr(t) — (Pt), are given in the frequency domain by 



8Pt(&) 



t^ePinSom^b ( ^) 

(rV2) 2 +A2 



(r - r e )/2 + iA (T -T e )/2- iA 



r,/2-i(Ao + fl) T t /2 + i(Ao-a) 



(D3) 



where 8Pt(Q) is the Fourier transform of 8Pr(t). By using Eq. (D3 ), we obtain a power spectral 
density (PSD) for the cavity transmission of 



s P (a) = g 2 om pls Xb (a)H(a), 



(D4) 
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where S Xb (£2) is the spectral intensity of the mechanical displacement for the optically bright mode 
which will be discussed in detail in the following sections. H(Q.) is the cavity transfer function 
defined as 

H( r« = T * 4A o( r o + ^ 2 ) m5 , 

[ [A 2 + (r f lif\ 2 + ^) 2 + (r t im [(A - a)2 + (r f /2)2] • 1 ; 

In general, when compared with S Xb (Q), H(Q.) is a slowly varying function of Q. and can be well 
approximated by its value at the mechanical resonance: H(Q.) H{£l' mb ). Clearly then, the power 
spectral density of the cavity transmission is linearly proportional to the spectral intensity of the 
mechanical displacement of the optically bright mode. 

APPENDIX E: THE MECHANICAL RESPONSE WITH MULTIPLE EXCITATION PATH- 
WAYS 

When the optically bright mode is coupled to an optically dark mode, the Hamiltonian for the 
coupled mechanical system is given by the general form: 

J%n=P- + + ^r~ + ^k d x 2 d + KX h X d , (El) 

where Xj, pj, kj, and mj (j = b, d) are the mechanical displacement, kinetic momentum, the spring 
constant, and the effective motional mass for the j th mechanical mode, respectively, and K repre- 
sents the mechanical coupling between the bright and dark modes. The subscripts b and d denote 
the optically bright and optically dark modes, respectively. With this system Hamiltonian, in- 
cluding the optical gradient force on the optically bright mode and counting in the mechanical 
dissipation induced by the thermal mechanical reservoir, we obtain the equations of motion for the 
two mechanical modes: 

d 2 Xb dxb 2 K Fb F a 

-nr + r mfc — + Q. mb xb H *d = 1 , (E2) 

dt z dt mi, m b mb 

d 2 x d dx d 2 K F d m ~ 

+ r md — + a md x d + -** = -, (E3) 

where Q^j = ^ is the mechanical frequency for the / h mode. Fj (j = b,d) represents the 
Langevin forces from the thermal reservoir actuating the Brownian motion, with the following 
statistical properties in the frequency domain: 

(Fi(a u )Fj (a v ) ) = ImiTmikBTdijlTtd (Q. u - Q v ), (E4) 
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where i,j = b,d, T is the temperature and kg is the Boltzmann constant. Fi(Q) is the Fourier 
transform of FAt). 



InEq. (E2\,F 



gomW 

too 



represents the optical gradient force. From the previous section, we 



find that it is given by 



Fo(t) 



Som 



a o\ 2 + a.Q8a(t) + ao8a*(t)~\ . 



(E5) 



The first term is a static term which only changes the equilibrium position of the mechanical 
motion. It can be removed simply by shifting the mechanical displacement to be centered at the 
new equilibrium position. Therefore, we neglect this term in the following discussion. The second 



and third terms provide the dynamic optomechanical coupling. From Eq. (C5), the gradient force 
is found to be given in the frequency domain by 



F (a) = f (a)x b (a) 



2g 2 om \a \ 2 A x h (Cl) 



A 2 -a 2 + (T t /2) 2 + iT t a 



mo [(A + n) 2 + (r t /2) 2 } [(Ao - a) 2 + (ry 2)2; 

which is linearly proportional to the mechanical displacement of the optically bright mode. 



, (E6) 



Equations fE2| ) and fE3] ) can be solved easily in the frequency domain, in which the two equa- 
tions become 



L d (Q,)x d + 



K _ 

— *d = 
m b 


Fb F 
m b nib 


(E7) 


K _ 

Xb = 

m d 


m d '' 


(E8) 



where Lj(fl) = Q.f n : — Or — iT m jQ. (j = b,d). Substituting Eq. (E6) into Eq. (E7), we find that 



Eq. (E7 ) can be written in the simple form, 



L b {Ll)x b ^ x d = — . 



(E9) 



where L&(£2) is now defined with a new mechanical frequency Of, and energy decay rate T' mb as 



L b (Q) = Q mb — d} — iT mb ci 

and the new Ql mS and Y' mS are given by 

2 ^Ll«o| 2 A 



urn 



(a' mb ) 2 -a 2 -ir' mb a, 



(E10) 



(PLY 



A 2 -n 2 + (T t /2) 2 



Q2mb+ m b (Oo [(A + a)2 + (r f /2)2][(A -a)2 + (r,/2)2] 



2 , 2gSml«o| Ao 
*> L mb ' 



^-^L + (r./2) 2 



^'mb = E 



mb 



m b COo [(Ao + Qmb) 2 + (T t /2) 2 ] [(Ao - a mb ) 2 + (T t /2) 2 ] 

2^2 m l«ol 2 r f A 1 



(Ell) 



m b coo [(A + a)2 + (r f /2)2] [(A - Q) 2 + (F t /2) 2 } 



^mb 



1 



m b (Oo [(A + tomb) 2 + (IV 2)2] [(A - a mb ) 2 + (iy 2)2] 
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(E12) 



Clearly, the effect of the optical gradient force on the optically bright mode is primarily to change 
its mechanical frequency (the optical spring effect) and energy decay rate (mechanical amplifica- 
tion or damping). 



Equations (E8) and (E9) can be solved easily to obtain the solution for the optically bright 
mode, 



x b (Ci) 



L d {Q) 



m b m d 



L b (a)L d (a)-^ 



(E13) 



where r\ L 



m b m d 



represents the mechanical coupling coefficient. By using Eq. (E4) and (E13), 



we obtain the spectral intensity of the mechanical displacement for the optically bright mode, 



S Xb (&) = 



2k B T r) A T md + T mb \L d (Q. ) 
m b \L b {n)L d {Q)-T\* 



2 ' 



(E14) 



where L b (Q.) is given by Eq. (E10). The mechanical response given by Eq. (E14) is very similar 
to the atomic response in EIT. 

As discussed in the previous section, the power spectral density (PSD) of the cavity trans- 
mission is linearly proportional to Eq. (E14). Equation ( E14| ) together with ( |D4[ ) is used to 
find the theoretical PSD shown in Fig. 3, by using an optomechanical coupling coefficient of 
g om /2n = 33 GHz/nm and an effective mass of m b = 264 pg for the flapping mode, both obtained 
from FEM simulations. The intrinsic and loaded optical quality factors of 1.07 x 10 6 and 0.7 x 10 6 
are obtained from optical characterization of the cavity resonance, and are also given in the cap- 
tion of Fig. 1 of the main text. The intrinsic mechanical frequencies and damping rates of the two 
modes (Q. mb , 0. md , T mb , and T md ) are obtained from the experimentally recorded PSD of cavity 
transmission with a large laser-cavity detuning, as given in the caption of Fig. 3 of the main text. 
The mechanical coupling coefficient r\ is treated as a fitting parameter. Fitting of the PSDs results 
in 7] = 3.32 MHz, indicating a strong internal coupling between the two mechanical modes. As 
shown clearly in Fig. 3d, f-h of the main text, our theory provides an excellent description of the 
observed phenomena. 



APPENDIX F: THE MECHANICAL RESPONSE WITH EXTERNAL OPTICAL EXCITATION 

The previous section focuses on the case in which the mechanical excitations are primarily in- 
troduced by the thermal perturbations from the environmental reservoir. However, the mechanical 
motion can be excited more intensely through the optical force by modulating the incident optical 
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wave. In this case, the input optical wave is composed of an intense CW beam together with a 



small modulation: A,-„ = A ; „o + 8A(t). As a result, Eq. (C3 ) now becomes 

dda 



(jA - T t /2)8a - ig om x b a Q + i^Y e 8A. 



dt 



(Fl) 



This equation leads to the intracavity field modulation given in the frequency domain as: 

igomaox b (Q.) - iy / T\8A(Q.) 



/(Ao + fl)-r,/2 



(F2) 



where 8A(£i) is the Fourier transform of 8A(t). By use of this solution together with Eq. (E5 ), the 
gradient force now becomes 

f (q) = f (a)x b (a)+F e (a), (F3) 



where f (Q) is given by Eq. (E6) and F e (Q) represents the force component introduced by the 



input modulation. It is given by the following form: 



i\/T~ P 



(Oq 



a 8A*(-Q) 
i(A + £2)-r,/2 ' i(A-Q.) + F t /2 



a* 8A(£l) 



+ 7 



(F4) 



In particular, in the sideband-unresolved regime, Eq. (F4) can be well approximated by 



Fe(&) 



a* SA(n)+a SA*(-£l) 



(F5) 



CBb(iA-r,/2) 

In the case that the mechanical excitation is dominated by the external optical modulation, the 
thermal excitation from the reservoir is negligible and Eqs. ( |E2[ ) and ( |E3[ ) become 



d 2 x b 



dx b 



K 



dt 2 + Tmb ^i + QmbXb + m Xd ~m' 
d 2 x d dx d 2 K n 

- w+ r md —+n m , Xd +-x b = o. 



(F6) 
(F7) 



Using Eqs. ( |F2| ) and ( [F3] ), following a similar procedure as the previous section, we find that the 
mechanical displacement for the optically bright mode is now given by 



(F8) 



m b L b (Cl)L d (Cl)-r)^ 

where L b {Q) and L d (Cl) are given in the previous section. Clearly, the mechanical response given 



in Eq. (F8 ) is directly analogous to the atomic response in EIT systems [38J. 
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APPENDIX G: CORRESPONDENCE OF CAVITY OPTOMECHANICS TO COHERENT 
STOKES AND ANTI- STOKES RAMAN SCATTERING 



In this section, we show a direct correspondence between cavity optomechanics and coherent 
Stokes and anti-Stokes Raman scattering. The system Hamiltonian of an optomechanical cavity is 
given by the following general form: 

M' = hojQa^a + hQ. m b^b + hg om Xba^a, (Gl) 

where a and b are the annihilation operators for photon and phonon, respectively, normalized such 
that a* a and b^b represent the operators for photon and phonon number. x\, is the mechanical 
displacement for the optically bright mode, related to b by 



x b 



n 



b+bA. (G2) 



Therefore, the interaction Hamiltonian between the optical wave and the mechanical motion is 
given by 

J% = hga*a(b + b^ , (G3) 

( s 2 Ti 3 \ 1 / 2 

where the factor g= [fi^ ) . 

The mechanical motion modulates the intracavity field to create two optical sidebands. As a 
result, the optical field can be written as 

a = a p + a s e~ i£imbt + aje i£1 "' h ' ', (G4) 

where a p is the field amplitude of the fundamental wave, and a s and a, are those of the generated 
Stokes and anti-Stokes wave, respectively. As the magnitudes of the Stokes and anti-Stokes side- 



bands are much smaller than the fundamental wave, when we substitute Eq. (G4) into Eq. (G3 ) and 
leave only the first-order terms of a s and a,;, under the rotating-wave approximation, the interaction 
Hamiltonian becomes 

= hg(b + b^ apcip + tigb** ^fljflp + aj/Z;- J + tigb {^ p a s + ci]a p ^ . (G5) 



In Eq. (G5), the first term describes the static mechanical actuation, which changes only the equi- 
librium position of mechanical motion and is neglected in the following analysis, as discussed 
previously. The second and third terms show clearly that the process corresponds directly to co- 
herent Stokes and anti-Stokes Raman scattering as shown in Fig. 4d in the main text. 
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APPENDIX H: THE MECHANICAL RESPONSE WITH THREE MODE COUPLING 



The coherent mixing of mechanical excitation is universal to gradient-force-based NOMS with 
a giant optical spring effect. Similar phenomena to that presented for double-disks were also 
observed in the coupled photonic crystal nanobeams. However, due to its device geometry, the 
coupled nanobeam has more complex mechanical mode families in which all the even-order me- 
chanical modes are optically dark, because they exhibit a mechanical node at the beam center 
where the optical mode is located. As the same-order common and differential motions of the two 
beams have similar mechanical frequencies, they can simultaneously couple to an optically bright 
mode, leading to multiple excitation interferences on the mechanical response. 

In the case when the optically bright mode is coupled to two optically dark modes, the Hamil- 
tonian for the mechanical system is given by the following general form: 

^ m= L (^ + ^k i xf\+K l X b Xi + K 2 X b X2 1 (HI) 
i=b,l,2 \ 1 ' 

where i = b, 1 , 2 corresponds to the optically bright mode and optically dark modes 1 and 2, re- 
spectively. With this Hamiltonian, counting in both the optical gradient force and the Langevin 
forces from the thermal reservoir, we obtain the equations of motions for the three modes: 



d 2 x b ^b, 2 , *i , K 2 Fb , F 
dt 2 + mb dt mb 



+ F mb — + & mh Xb H xi H x 2 = 1 , (H2) 

m b m b m b m b 



d 2 X{ dx\ 2 K i F\ 

- 1 - 2 -+r ml — + Ll ml x 2 ^ x b = — , (H3) 

dt z dt mi mi 



d 2 x 2 1—, dx 2 /~\2 K 2 F 2 



mi mi 

+ r w2 =f + ^ 2 2 x 3 + —*b = — , (H4) 
dt L dt m 2 m 2 

where the gradient force F a is given by Eq. ( |E5[ ), and the statistical properties of the Langevin 



forces are given by Eq. (E4). Following the same analysis as the previous section, we can obtain 



the spectral intensity for the mechanical displacement of the optically bright mode as 

2^rT ? 4r ml |L 2 (a)| 2 + 77 2 4 r m2 |L 1 (^)| 2 +r ra ,|L 1 (^)L 2 (^)| 2 

Sx b\ il ) = — r — 7Z — 7— ~2 ' ( H5 ^ 



2k B T 7] A x Y m 


l\Li{a)\ 2 + l\iT„a 


L l (o.)\ 2 +r mb 


L!(^)L 2 (a)| 2 


m b 


L b {Q)Li{Q)Li{Q) 




2 



K 2 

where rjj = -^j^- (j = 1,2) represents the mechanical coupling coefficient. Lj(Cl) = £l 2 nj , — £l 2 



1 J m b mj yJ ' ' 1 ~r «= 7V~y mj 



ir m jQ, (j =1,2) and L b (Q.) is given by Eq. (E10) with Q,' . and T' mb given in Eqs. (Ell) and 



( |E12 ), respectively. As the optical wave is coupled to the optically bright mode only, the power 



spectral density of the cavity transmission is still given by Eq. (D4), with the mechanical response 



S Xb given in Eq. (H5) 
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FIG. 5; (a) Experimentally recorded power spectral densities of the cavity transmission for the zipper cavity described in the main text, with 
an input power of 5.1 mW. Each curve corresponds to a laser frequency detuning indicated in (c). Each curve is relatively shifted by 5 dB in 
the vertical axis for a better vision of the mechanical frequency tuning and the induced mechanical interference. The optically dark mode II and 
III have a full-width at half maximum (FWHM) of 0.16 and 0.15 MHz, respectively. The optically bright mode I has an intrinsic FWHM of 
0.30 MHz. (b) The corresponding theoretical spectra of the power spectral density, (d)-(f) The detailed spectra of the power spectral density at 
three frequency detunings indicated by the three arrows in (c). The blue and red curves show the experimental and theoretical spectra, respectively, 
(h) FEM simulated mechanical motions for the fundamental differential mode (I), the second-order common (II) and differential (III) modes, whose 
frequencies are indicated by the arrows in (a). The color map indicates the relative magnitude (exaggerated) of the mechanical displacement. 



Figure [5] shows the PSD of the cavity transmission by launching a continuous wave into a reso- 
nance of the coupled nanobeams with an intrinsic and loaded Q factor of 3.0 x 10 4 and 2.8 x 10 4 , 
respectively. Three mechanical modes are clearly visible, where mode I is the fundamental dif- 
ferential mode [Fig.[5pi)I], and mode II and III correspond to the second-order common and dif- 
ferential modes [Fig.(5fh)II and III], respectively. Similar to the double-disk NOMS, the gigantic 
optical spring effect shifts the frequency of the optically bright mode I from its intrinsic value of 
8.06 MHz to 19 MHz, crossing over both optically dark modes II and III closely located at 16.54 
and 17.04 MHz and resulting in complex interferences on the power spectra [Fig. |5ja)]. Equation 
(H5) provides an accurate description of the observed phenomena, as shown clearly in Fig.^b), 
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(d)-(f). Fitting of the PSD results in mechanical coupling coefficients of 77 1 =3.45 MHz and 
rj2 = 3.48 MHz, implying that the two optically dark modes couple to the fundamental optically 
bright mode with a similar magnitude. 



APPENDIX I: RENORMALIZATION OF MECHANICAL MODES BY THE GRADIENT 
FORCE 

In addition to the gradient force and Langevin force, mechanical motion of the two microdisks 
(or two nanobeams) changes the gap between them and thus introduces a pressure differential be- 
tween the gap and the outer region [[32l . which functions as a viscous force to damp the mechanical 
motion. In the previous sections, we have incorporated this squeeze-film effect in the mechanical 
damping rate. However, as the generated pressure differential is sensitive to the gap variations, 
we can expect that the squeeze-film effect behaves quite differently when the two disks/beams 
vibrate independently or cooperatively. Therefore, we describe it explicitly here for the analysis 
of mechanical mode renormalization. 

In general, the motion of individual disks or nanobeams satisfies the following equations: 

d 2 xi dxi 2 F x F F q 

+r m i— +£2 ml xi = — 1 1 — , (ii) 



mi m\ m\ 



dt 2 ' mJ dt 
d 2 x 2 dx 2 2 F 2 F F q 

-Hi- + T 'n2 —J— + ^ m 2*2 = " , (12) 

dt z dt tn.2 ni2 ni2 

where F q is the viscous force from the squeeze film damping, and rrij, Xj, Cl m j, T m j, Fj (j = 1,2) 
are the effective mass, the mechanical displacement, resonance frequency, damping rate, and the 
Langevin force for individual disks (or beams), respectively. 

The optically bright mechanical mode corresponds to the differential motion of the two 
disks/beams, with a mechanical displacement given by Xb = x\ — X2- By transferring Eqs. ( pTj ) 



and (12) into the frequency domain, it is easy to find that the mechanical displacement of the 



optically bright mode is given by 

m\L\{Ll) m.2L2(il) 



1 1 

+ ■ 



m\L\(Q) m,2L2(£l) 



FJQ)+F (£l) 



(13) 



where Lj(Q.) = Q. 2 mj - Q 2 - iT mj Q. (j = 1,2). The squeeze-film effect is produced by the pressure 
differential between the gap and the outer region introduced by the differential mechanical motion, 
and thus has a magnitude linearly proportional to the differential displacement. In general, it can 
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be described by F q (Ci) = f q (Q.)xb(Q) , where f q (£l) represents the spectral response of the squeeze 
gas film [32]. Using this form together with Eq. fEo} in Eq. ([D]), we obtain the spectral intensity 
of the optically bright mode displacement, 



S*(0) 



2k B T 



Sf|L 2 (a) | 2 + ^|Li(n)| 



Li(n)L2(Q)-[/ (Q)+/ 9 (a 



Li(fl) L 2 (H) 



2' 



(14) 



m2 mi 

As the squeeze-film effect primarily damps the differential motion, its spectral response can 
be approximated as f q (Q.) ~ ia 9 G. Moreover, since the two disks or nanobeams generally have 
only slight asymmetry due to fabrication imperfections, they generally have quite close effective 
masses and energy damping rates: m\ m m 2 = 2m b and T m \ « T m 2 = T m , where we have used the 
fact that the effective motional mass of the differential motion is given by nib = mimij (m\ + m%) . 



As a result, Eq. (|I4J) can be well approximated by 



Li(n)| 2 + |ia(Q)| 2 



m b |Li(ft)L2(fl)-i [/ (a)/m fc + iT 9 ft] [Li(fi)+I^(Q)]| 



(15) 



where T q = OC q /mb represents the damping rate introduced by the squeeze gas film, and the spectral 



response of the gradient force f (Cl) is given by Eq. (E6 ). 

The intrinsic mechanical frequencies of 7.790 and 7.995 MHz for the two individual nanobeams 
are measured from the experimental recorded PSD with a large laser-cavity detuning. The optome- 
chanical coupling coefficient is 68 GHz/nm and the effective mass is 10.75 pg for the fundamental 
differential mode, both obtained from FEM simulations (note that these values are different than 
those quoted in Ref. 11221 due to the different definition of mode amplitude for Xb). The intrinsic 
and loaded optical Q factors are 3.0 x 10 4 and 2.8 x 10 4 , respectively, obtained from optical char- 
acterization of the cavity resonance. By using these values in Eqs. ( [15] ) and (E6), we can easily 
find the mechanical frequencies and linewidths for the two renormalized modes, where we treat 
the intrinsic mechanical damping rate Y m and the squeeze-film-induced damping rate Y q as fitting 
parameters. As shown in Fig. 2 of the main text, the theory provides an accurate description of the 
mechanical mode renormalization, with a fitted intrinsic mechanical and squeeze-film damping 
rate of 0.03 and 0.2 MHz, respectively. 

Similarly, we can obtain the spectral intensity of xj =x\+X2 for the optically-dark mechanical 
mode, which is given by the following form: 



S x iQ)=2k B T 
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i*(G)-£[/ (G)+/,(n 



L 1 (Q)L 2 (G)-[/ (n)+/ 9 (Q)] 



Li(Q) , Zq(Q) 
ni2 mi 



(16) 
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Similar to the optically-bright mode, with m\ « ni2 = 2m b and r m i « T ra 2 = r„„ Eq. ([16]) can be 
well approximated by 

Sx (a) « fcg7Tm l^i(^)-M^)l 2 + l^(a)-M^)l 2 ^ (I7) 

m ^ |Li(^)L2(^)-^(a)[Li(^)+L 2 (^)]| 2 ' 

where = m/2 is the effective mass of the common mode and h(Q.) = [/ (Q)/m& + jT 9 Q] 
represents the total spectral response of the optical gradient force and squeeze film damping. In 
particular, when the optical-spring-induced frequency shift is much larger than the intrinsic me- 
chanical frequency splitting, the spectral intensities of these two modes reduce to 



MQ)« ^TT m /m b 2k B TT m /m d 

\L (Q)-h(a)\ 2 May 2 



where Lq(CI) = {0. m \ + £l m 2) /4 — Or — iT m Q.. Equation (18) indicates that the optically bright 
and dark modes reduce to a pure differential and common modes, respectively. 
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